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Abstract 


\ 

A  most  favorably  oriented  crystal  located  at  a  free  surface  of  a  f.c.c.  polyciystal  sub¬ 
ject  to  creep  under  cyclic  tension  and  compression  of  high-cycle  fatigue  is  considered.  An 
extrusion  in  this  crystal  is  shown  to  be  produced  by  a  positive  slip  in  one  thin  slice  "P"  and  a 

r tnj. 

negative  slip  in  a  closely  located  slice  ”Q".  An  initial  tensile  strain  s^in  the  thin  slice  R 
sandwiched  between  P  and  Q  causes  a  positive  initial  shear  stress  T^in  P  and  a  negative  one 
in  Q.  It  is  shown  that  the  extrusion  growth  causes  a  tensile  strain  in  R,  which  can  activate  a 
second  slip  system  giving  a  creep  strain  with  a  tensor  component's,^  It  has  the  same  effect 
as  the  initial  strain  'e^in  causing  this  difference  in  shear  stresses  in  P  and  Q  and  gives  much 
additional  extrusion  growth.  The  extent  of  intrusion  and  extrusion  is  important  in  this  study 


of  crack  initiation. 
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Introduction 


Single  crystal  test  at  room  temperature  (Taylor,  1938)  have  shown  that,  under  stress, 
slip  occurs  along  certain  directions  on  certain  planes.  Slip,  the  result  of  displacement  of  dislo¬ 
cations.  has  been  found  to  depend  on  the  resolved  shear  stress  and  not  on  the  normal  stress  on 
the  sliding  plane.  This  dependency  of  slip  on  the  resolved  shear  stress,  known  as  Schmid’s 
law,  has  been  shown  by  Parker,  1961  to  hold  also  under  cyclic  loadings.  Deformation  of  sin¬ 
gle  crystals  at  elevated  temperatures  (Johnson,  et  al,  1953,  1955)  also  occurs  by  slip  in  pri¬ 
marily  the  same  slip  systems  that  are  operative  at  room  temperature.  This  slip  rate  was  found 
to  depend  on  the  resolved  shear  stress.  This  agrees  with  Gilman’s  1969  observation  that  this 
dislocation  velocity  varies  with  the  resolved  shear  stress. 

Following  the  clue  provided  by  the  observation  of  extrusions  and  intrusions  in  fatigue 
specimens  (Forsyth,  1954),  a  number  of  distinguished  investigators;  Mott  1958,  Cottrell  and 
Hull  1957,  Thompson  1959,  McEviley  and  Machlin  1959,  Wood  1956  and  others  have  pro¬ 
posed  different  theories  of  fatigue  crack  initiation.  The  above  early  fatigue  crack  initiation 
theories  primarily  consider  the  paths  of  dislocation  movement  without  considering  the 
resolved  shear  stress  induced  by  this  movement.  This  seems  to  be  the  main  cause  for  these 
theories  to  be  not  quite  satisfactory,  as  pointed  out  by  Kennedy  1961.  Lin  and  Ito  1969  pro¬ 
posed  first  quantitative  micromechanic  theory  of  fatigue  crack  initiation.  This  theory  calcu¬ 
lates  the  stress  fields  produced  by  slip  in  fatigue  bands  and  has  been  shown  to  agree  with 
nearly  all  experimental  observations  (Lin,  1977,  1981).  This  theory  is  briefly  reviewed  here 
and  is  applied  to  the  present  analysis. 
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Initial  defects  always  exist  in  metals  and  cause  an  initial  stress  field.  Let  the  initial 
resolved  shear  stress  be  denoted  by  T7.  During  loading,  when  the  resolved  shear  stress  in 
some  region  reaches  the  critical  shear  stress  zc,  slip  occurs.  After  unloading,  this  slip  remains 
and  induces  a  residual  resolved  shear  stress  t*.  The  resolved  shear  stress  due  to  loading  is 
denoting  by  T'4.  The  total  shear  stress  after  reloading  is 


t  =  X7  +  X*  4-  T4 


(1) 


Consider  a  most  favorably  oriented  crystal  at  a  free  surface  of  a  f.c.c.  polycrystal  as 
shown  in  Fig.  1.  For  an  extrusion  to  start  in  a  thin  slice  R  sandwiched  between  two  slices  P 
and  Q,  positive  shear  must  occur  in  P  and  negative  in  Q.  The  initial  shear  stress  field  x7, 
favorable  for  this  sequence  of  slip,  clearly  is  one  having  positive  shear  stress  in  P  and  nega¬ 
tive  in  Q.  Such  an  initial  stress  field  can  be  provided  by  an  initial  tensile  strain  e^a  *n  R- 
This  positive  e4a  can  be  provided  by  a  row  of  interstitial  dislocation  dipoles  and  a  negative 
£aa  by  vacancy  dislocation  dipoles  as  shown  by  Lin  and  Ito  1969. 


A  tensile  loading  x22  on  the  polycrystal  (Fig.  1)  produces  a  positive  T4  in  the  whole 
crystal.  Taking  x7  to  be  positive  in  P  and  negative  in  Q,  we  have  x7  +  X4  in  P  reaching  the 
critical  shear  stress  T?  first;  and  hence,  P  slides.  Due  to  the  continuity  of  stress  field,  slip  in  P 
relieves  not  only  the  positive  shear  stress  in  P  but  also  in  Q  (Lin  and  Ito,  1969,  Lin  1977). 
Hence,  this  slip  increases  the  negative  resolved  shear  stress  in  Q  to  cause  Q  to  slide  more 
readily  in  the  reverse  loading.  The  negative  slip  in  Q  relieves  the  negative  shear  stress  not 
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Figure  1.  Most  Favorably  Oriented 
Crystal  at  Free  surface 


only  in  Q,  but  also  in  P,  thus  causing  P  more  readily  to  slide  in  the  next  forward  loading. 
This  process  is  repeated  for  every  cycle  thus  providing  a  natural  gating  mechanism  for  a 
monotonic  buildup  of  local  slip  strain  e^p  in  P  and  Q,  pushing  R  out  of  the  free  surface  and 
starts  an  extrusion.  Interchanging  the  signs  of  the  initial  stresses  in  P  and  Q  initiates  an  intru¬ 
sion  instead  of  an  extrusion.  This  theory  is  extensively  supported  by  metallurgical  observa¬ 
tions  (Lin,  1977). 

Extent  of  Extrusion 

The  buildup  of  this  local  slip  strain  e^p  in  P  and  Q  is  caused  by  the  positive  and  nega¬ 
tive  initial  shear  stresses  x^p  which,  in  turn,  is  caused  by  in  R.  If  R  were  cut  out,  this 
free  length  of  R  would  be  longer  than  the  slot  cut  by  an  amount  referred  to  as  the  "static 
extrusion"  (Mughrabi  et  al,  1983).  This  e/aa  causes  an  initial  compression  x^a  in  R.  As  the 
extrusion  grows  under  cyclic  loading,  the  slice  R  increases  in  length.  This  elongation  causes 
the  compression  to  decrease.  A  question  has  been  raised  as  to  whether  the  extrusion  growth 
will  cease  after  the  extrusion  has  reached  the  static  extrusion.  Lin,  et  al.,  1988  have  shown 
that  the  residual  tensile  stress  x£a  caused  by  elongation  in  R  due  to  extrusion  can  cause 
changes  of  resolved  shear  stresses  in  all  twelve  slip  systems.  The  resolved  shear  stress  in  one 
slip  system  may  reach  the  critical  and  slide.  The  creep  strain  caused  by  slip  in  this  sys¬ 
tem  has  a  tensor  component  e£a  just  like  e^a  in  causing  the  positive  and  negative  x^p,  respec¬ 
tively,  in  P  and  Q.  This  secondary  slip  has  been  shown  by  Lin,  et  al.,  1988  to  increase 
greatly  the  extent  of  extrusion  and  intrusion  in  time-independent  slip.  This  paper  shows  the 
this  secondary  slip  also  increases  this  extent  of  intrusion  and  extrusion  in  fatigue  with  creep. 
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Residual  Stress  Field  Caused  by  Creep  Strain 

The  creep  rate  depends  on  the  resolved  shear  stress  x  which  is  given  in  Eq.  1.  as  the 
sum  of  the  initial,  applied  and  residual  shear  stresses.  To  calculate  x,  the  analogy  of  creep 
strain  and  applied  force  is  applied.  It  has  been  shown  that  the  equivalent  body  force  per  unit 
volume  along  xl  -  axis  due  to  creep  strain  is 

F,  =  -Q&ij  £ckjc,i  +  IGtfjj) 

~  (2) 

where  X  and  G  are  Lame’s  constants  and  8^  is  the  Kronecker  delta.  The  repetition  of  sub¬ 
script  denotes  summation  and  the  subscript  after  comma  denotes  differentiation,  with  respect 
to  the  coordinate  variable.  The  equivalent  surface  force  per  unit  area  along  the  x,  -  axis  also 
has  been  shown  as 

Si  =  (X£^5,  +  2 G  ef  )  v- 

~  (3) 

where  vy  is  the  cosine  of  the  angle  between  the  exterior  normal  to  the  surface  and  the 
Xj  -  axis. 

From  a  given  e-jj,  we  obtain  F,  and  5,  from  Eqs.  (2)  and  (3).  Then  we  apply  this  F,- 
and  5,  to  the  body  and  obtain  the  strain  distribution  e,y  in  this  body  by  elastic  solution.  The 
residual  stress  field  in  this  body  is  then 

xfj  =  Xdtj  e**  +  2|i£,y  -  X5y  -  2[L  efy 

(4) 

This  analogy  becomes  Duhamel’s  analogy  for  thermal  stress  (Lin,  1968),  if  the  creep  strain  efj 
is  replaced  by  thermal  strain  ejj. 
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The  lengths  of  the  strips  P,  Q  and  R  along  the  j^-direction  are  much  larger  than  the 
thickness:  the  thickness  also  is  much  smaller  than  their  inclined  lengths  shown  in  Fig.  1. 
Hence,  strains  in  the  major  central  portion  of  the  strips  are  taken  to  be  of  generalized  plane 


strain,  i.e. 


M;  =  *2)  i  =  1,2,3 


This  gives 


xij  2w  y_2u  +  2  + 


where  0  =  Mltl  +  m2i2  and  v  is  Poisson’s  ration.  Substitution  of  the  expression  into  the  condi¬ 
tions  of  equilibrium  yields 

V  ua  +  TT-  3 —  +  -77*  =  0*  a  =1,2 

l-2u  dxa  G  m 


where 


V2m,  +  —  =0 


v2  =  JL  +  JL 

dxj  d*2 


where  Fa  and  F3  can  be  either  the  applied  forces  or  the  equivalent  forces. 


The  above  two  equations  (7)  and  (8)  are  not  coupled  and  can  be  solved  separately. 


Lin  and  Lin  1974  and  Lin,  et  al.,  1988  have  shown  the  solution  of  this  stress  field  caused  by 
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Fa  by  using  Airy’s  stress  functions.  The  solution,  of  the  displacement  fields  and  then  the 
stress  fields  caused  by  F3  also  has  been  shown  by  Lin,  et  al.,  198R.  Consequently  the  stress 
fields  caused  by  the  equivalent  forces  due  to  creep  strain  can  readily  be  calculated. 

The  twelve  slip  systems  of  a  f.c.c.  crystal  are  shown  in  Fig.  2.  The  most  favorable 
slip  system  in  the  most  favorably  oriented  crystal  making  45°  with  the  direction  of  loading  is 
identified  with  (Fig.  2).  During  fatigue  loading,  this  buildup  of  large  local  creep  shear 
strain  eap  in  P  and  Q  causes  R  to  elongate,  inducing  an  appreciable  direct  stress  It  has 
been  shown  by  Lin,  et  al.,  (1988)  that  there  are  four  slip  systems  c}i  c3  dl  and  d2  equally 
favorable  under  x^.  Among  the  four  c3  is  the  most  favorable  under  this  applied  load  x22. 
Hence,  c3  is  considered  to  be  active  in  R. 


Figure  2.  Crystallographic  Directions  of  a  f.c.c.  Crystal 
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The  normal  to  the  slip  plane  P  and  the  slip  direction  a  of  the  slip  system  a2  lie  in  the 
as  plane.  The  equivalent  forces  Fa  caused  by  slip  in  the  a2  -  slip  system  are  in  this 
plane.  Hence  plane  strain  solution  satisfying  the  differential  equation  (7)  was  used  (Lin  and 
Lin  1974).  However  the  slip  direction  q  and  the  normal  to  the  slip  plane  T)  of  the  c3  slip  sys¬ 
tem  are  not  in  the  xx  x2  -  plane.  The  equivalent  forces  caused  by  slip  in  the  c3  slip  sys¬ 
tem  does  not  lie  in  the  x2  -  plane.  Hence  the  equivalent  force  has  a  component  f3.  Then 
the  stress  fields  caused  by  the  components  F3  lying  in  the  plane  are  solved  by  Eq.  7  and  those 
caused  by  the  component  F2  are  solved  by  Eq.  8.  (Lin,  et  ah,  1988).  This  shows  how  this 
residual  stress  caused  by  creep  in  the  primary  slip  system  e„n  and  those  in  the  secondary  slip 
system  are  analytically  calculated. 

Numerical  Calculations 

The  dimensions  of  the  most  favorably  oriented  crystal  at  the  free  surface  of  a  polycry¬ 
stal  of  pure  aluminum  are  shown  in  Fig.  1.  The  thickness  in  P  and  Q  is  0.01  |i,  and  the 
thickness  in  R  is  0.  ’  |i  along  the  x2  direction.  As  shown  in  Fig.  4,  the  initial  shear  stress  in 
P,  Q  is  taken  to  vary  linearly  from  0  at  the  interior  boundary  to  1.1  x  104  Pa  (1.6  psi)  at  the 
free  surface.  This  is  expressed  as 

4(a)  =  1.1  x  10*(1  - 

From  equilibrium  conditions 

a  2 
C(o)  -  -2 KuW*  =  -2.2  X  104(a-  looax^) 
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£/fia(a)  ~  E  ]  Taa  "  °'412  X  l°  5(Ct  100  x  V2  1 

The  static  extrusion  is  then 
50V2 

J  ec«x(a)^a  =  0.687  x  10-2 

O 

The  creep  strain  at  this  free  surface  corresponding  to  the  above  static  extrusion  is  0.486. 
The  single  crystal  creep  property  is  taken  from  tests  at  204.4°C  (400°F)  performed  by  John¬ 
son  et  al,  1953,  1955  and  is  represented  approximately  by  a  linear  relation  (Lin  and  Lin,  1978, 
1980). 

e„p  =  3.26  x  10-7[-Caf}-1.3167  x  106A7m2(  190.1  /wi)]/min 

where  the  dot  on  top  denotes  the  time  derivative  and  Tap  is  the  resolved  shear  stress  in  N/m2 
in  the  active  slip  system.  Creep  strain  is  highly  localized  in  slip  bands.  The  creep  rate  in  this 
band  is  much  higher  than  the  average  rate  of  the  whole  crystal  and  is  assumed  to  be  1000 
times  that  of  the  crystal  giving 

=  3.26  x  10-4(xaa-1.3167  x  l06A7m2)/min 

(9) 

Referring  to  Fig.l,  the  thin  strips  P,  Q  and  R  are  divided  along  their  lengths  into  thin 
parallelogram  grids.  For  numerical  calculation,  the  creep  strain  in  each  grid  is  assumed  to  be 
constant.  Hence,  there  is  no  creep  strain  gradient  and  no  equivalent  body  force  in  each  grid. 
However  each  grid  is  subject  to  a  set  of  equivalent  surface  forces  on  the  boundary.  The  aver¬ 
age  residual  resolved  shear  stress  in  the  ith  grid  caused  by  constant  creep  strain  in  the  jth 
grid  or  in  the  kth  grid  was  then  calculated  using  the  generalized  plane  strain  semi-infinite 
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solid  solution  given  by  Eqs.  7  and  8  and  is  denoted  by  -C(i,  ap;  j,  aP)  or  -C(i,  a(3;  k,  qrO- 
With  e£p  occurring  only  in  P,  Q  and  in  R,  the  residual  resolved  shear  stress  in  the  ith  grid 
due  to  creep  strain  e^p  and  is  then  written  as 


=  -CO,  ap;  j,  aP)e«pJ  ~  C(i,  ap;  k,  ^rt)e|n* 
tfr,,  =  -C(i,  ^T|;  j,  ap)eU.  -  C(i,  ^r\;  k, 

(10) 

The  repetition  of  subscript  "j"  denotes  summation  covering  all  grids  with  e^p  and  "k"  summa¬ 
tion  covering  all  grids  with  e^.  The  resolved  shear  stress  in  the  ith  grid  is  written  as 

TaP,  =  ^p,  +  ^ap,  +  ^p,- 

(11) 

x$n.  =  +  +  Tfn,- 

From  Eqs.  9  through  1 1  we  can  write 

e^p.  =  3.26  x  lO^^p  -  1.3167  x  l06A7m2(  190.1  psi)/min 

=  3.26  x  lO^tx^p.  +  x^p.  -  C(i,  aP;  j,  aP)e£p 

-  C(i,  aP;  k,  %T[)Ztr\  -  1-3167  x  106A7m2(  190.1  psi)} 

(12) 

A  similar  equation  can  be  written  for  e^.. 

From  the  above,  the  creep  strain  rates  in  grids  with  resolved  shear  stress  greater  than 
this  critical  1.3167  x  106A7m2  (  190.1  psi)  can  be  determined. 
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Write  Eq.  12  in  matrix  form 


{£}  =  [A]{e}  + 

(13) 

subject  to  the  initial  condition  (e(r0)}  =  {£„}.  where 

Bt  =  3.26  x  lCT4^.  +  t£P;  -  1.316  x  106A7m2) 

(14) 

The  series  solution  of  Eq.  13  is  given  by  Frazer  et  al,  1963.  The  complementary  solu¬ 
tion  of  Eq.  13  is  given  as 

{Ec}  =  e[/41(,-°{eo} 


=  {[/]  +  [A](t  -ta)+  ■■■  +  -\[A]P(t  -tQ)P+  •  •  ■  ]}  {£j 

P ■  (15) 

where  £0  is  this  creep  strain  at  ta,  e,  the  base  of  the  natural  logarithm,  and  |T|  is  the  identify 

matrix.  The  particular  solution  is  obtained  by  multiplying  Eq.  13  by  e~ M*.  This  gives 

±{e-[A]t{Zp})  =  e-[A],{B{X) } 

and 

t 

{£p}  =  e+W’fe-W*  '  ‘\B(s)}ds 

to 

The  cyclic  loading  is  taken  to  be  rectangular:  i.e.,  tension  and  compression  are  applied  instan¬ 
taneously  and  remain  constant  during  each  half  cycle. 

{ep}  =-{[/]- 
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{ep}  =  {/(*  -  ta)  +  -±-[A]it  -  g2  +  ....  +  -L[A]P-\t-  tQf 

L\  p\ 


+  •  • •  }{B} 


(16) 


To  calculate  the  creep  strain  distribution  in  a  given  time  interval  0  to  T.  This  interval 
is  divided  into  N  steps  with  A t  =  TIN.  The  At  is  chosen  tc  be  small  enough,  so  the  maximum 
eigenvalue  of  [A]Ar  is  less  than  unity.  Then  the  matrix  [A]n  (r  -  ta)n  approaches  a  zero  matrix 
as  "n"  approaches  infinity.  The  series  solution  of  Eqs.  15  and  16  will  converge. 


For  each  time  step  Ar,  take  the  first  p-terms  of  {ec}  and  {ep}  to  obtain  the  creep 
strains.  The  smaller  the  Ar  and  the  more  terms  p  used,  the  more  accurate  will  be  the  results. 
In  the  previous  work  (Lin  and  Lin,  1979,  1982),  the  creep  strain  rate  was  assumed  to  be  con¬ 
stant  in  At.  The  solution  corresponds  to  p  =  1.  It  has  been  found  that  using  this  series  solu¬ 
tion  reduces  the  numerical  calculation  by  more  than  six  times  for  the  same  accuracy  in 
evaluating  the  creep  strain  distributions.  In  the  present  example,  the  magnitude  of  applied 
load  was  2.7326  x  106  Mmm2  (380psi)  giving  f4  =  1.366  x  106  N/m2  (190psi)  Two  loading 
frequencies  were  calculated;  one  is  0.345  cycles/min  and  the  other  is  345  cycle/min.  In  the 
high  frequency  loading  the  half  period  was  equally  divided  into  ten  time  increments  and  three 
terms  in  the  series  solution  were  used.  The  creep  strain  in  the  second  slip  system  in  R 
versus  the  distance  from  the  free  surface  at  different  cycles  of  loading  under  a  frequency  of 
345  cycles  per  minute  was  calculated  and  is  shown  in  Fig.  3.  The  creep  strain  in  the 
second  slip  system  causes  some  £„a  which  in  turn  cause  additional  positive  +T^jj  in  P  and 
negative  in  Q.  The  creep  rate  in  P  and  Q  depends  on  the  sum  of  the  initial  and  the  resi- 
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dual  shear  stresses  t^p  +  t£p  in  P  and  Q.  This  sum  at  different  stages  of  loading  is  shown  in 
Fig.  4.  Note  that  this  sum  decreases  considerably  slower  in  the  cases  with  the  secondary  slip 
£*„  than  the  case  without.  When  this  sum  approaches  zero,  the  extrusion  approaches  static 
extrusion  and  ceases  to  grow.  The  creep  shear  strain  distributions  e£p  in  P  and  Q  at  different 
cycles  of  loading  were  calculated  for  the  cases  with  and  without  secondary  slip  and  are  shown 
in  Fig.  5.  The  e^p’s  at  the  free  surface  versus  cycles  of  loading  n  are  shown  in  Fig.  6  this 
strain  versus  loading  time  N/f  is  shown  in  Fig.  7.  Note  that  with  no  secondary  slip,  the  extru¬ 
sion  growth  ceases  when  the  extrusion  approaches  the  static  extrusion.  The  creep  strain  e„p 
in  P  and  Q  at  the  free  surface,  which  represents  the  amount  of  extrusion  or  intrusion  with  the 
secondary  slip  is  much  higher  than  that  without  the  secondary  slip.  Calculations  also  have 
been  made  at  a  frequency  of  loading  of  10  cycles/min,  and  the  creep  strain  in  P  and  Q  versus 
N/f  have  been  found  to  be  essentially  the  same  as  that  at  345  cycles/min.  In  the  high  cycle 
range,  the  creep  strain  depends  mainly  on  the  time  of  loading.  Similar  calculations  were 
made  on  a  low  frequency  loading  0.345  cycles/min.  The  creep  strain  at  this  low  frequency 
depends  on  the  number  of  cycles  as  well  as  N/f  the  time  of  loading. 
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Conclusions 


The  amount  of  creep  shear  strains  e^p  in  P  and  Q  at  the  free  surface  is  a  measure  of 
the  amount  of  extrusion  and  intrusion  and,  hence,  is  taken  as  a  measure  of  fatigue  crack  initia¬ 
tion.  The  calculated  results  indicate  that  this  surface  creep  strain  depends  only  on  the  time  of 
loading  N/f  in  high  frequency  range  and  becomes  dependent  also  on  the  number  of  loading 
cycles  in  the  low  frequency  range. 

Creep  strain  in  the  secondary  slip  system  has  a  tensor  component  e^a  which  has  the 
same  effect  in  causing  t^p  in  P  and  Q  as  ^aa.  Hence  this  creep  strain  e^p  in  P  and  Q  with 
this  secondary  slip  is  much  larger  than  that  without.  Without  the  secondary  slip,  the  extrusion 
ceases  when  it  reaches  the  static  extrusion.  However,  with  the  secondary  slip,  extrusion 
grows  considerably  beyond  the  static  extrusion.  This  may  be  a  main  cause  for  the  extent  of 
extrusion  to  be  much  more  in  a  face-centered  polycrystal  than  that  in  a  hexagonal  polycrystal. 
The  stress  intensity  factor  of  crack  depends  on  the  extent  of  extrusion  on  intrusion.  Hence, 
this  extent  is  important  for  the  study  of  fatigue  crack  initiation. 


15 


Acknowledgments 


The  support  of  the  U.S.  Office  of  Naval  Research  through  Contract  N00014-86-K-0153 
and  the  interest  of  the  scientific  officer.  Dr.  Yapa  Rajapakse,  are  gratefully  acknowledge. 


16 


References 


1.  Cottrell,  A.H.  and  D.  Hull,  "Extrusion  and  Intrusions  by  Cyclic  slip  in  Copper",  Proc. 
Roy.  Soc.  (London),  Vol.  242A,  p.  211,  1957. 

2.  Forsyth,  P.J.E.  and  C.A.  Stubbington,  "The  Slip  Band  Extrusion  Effect  Observed  in 
Some  Aluminum  Alloys  Subjected  to  Cyclic  Stresses",  J.  Inst.  Metals,  Vol.  83,  p.  395, 
1955. 

3.  Frazer,  R.A.  W.J.  Duncan  and  A.R.  Collar,  "Elementary  Matrices  and  Some  Applica¬ 
tions  to  Dynamics  and  Differential  Equations,"  Cambridge  University  Press,  p.  45, 
1938. 

4.  Gilman,  J.J.,  "Micromechanics  of  Flow  in  Solids",  McGraw-Hill,  New  York,  pp.  176- 
179,  1969. 

5.  Johnson,  R,D.,  R.R.  Shober  and  A.D.  Schope,  "The  Creep  of  Single  Crystals  of  Alumi¬ 
num",  National  Advisory  Committee  for  Aeronautics,  Technical  Note  2945,  1953. 

6.  Johnson,  R,D.,  A.P.  Young  and  A.D.  Schope,  "Plastic  Deformation  of  Aluminum  Sin¬ 
gle  Crystals  at  Elevated  Temperatures”,  National  Advisory  Committee  for  Aeronautics, 
Technical  Note  3351,  1955. 

7.  Kennedy,  A.J.  "Process  of  creep  and  Fatigue  of  Metals"  John  Wiley  &  Sons  Inc.  pp. 
331-343,  1963. 

8.  Lin.  T.H.,  "Micromechanic  of  Deformation  of  Slip  Bands  Under  Monotonic  and  Cyclic 
Loadings",  Reviews  of  the  Deformation  Behavior  of  Materials,  P.  Felham  (ed.),  Freund 
Publishing  House;  Tel-Aviv,  Israel,  1977. 

9.  Lin.  T.H.,  "Micromechanic  of  Fatigue  Crack  Initiation:  Theory  and  Experimental 
Observation",  Mechanics  of  Fatigue  -  AMD,  Vol.  47,  ASME,  Edited  by  T.  Mura,  pp. 
41-109,  1981. 

10.  Lin,  T.H.  and  Y.M.  Ito,  "Mechanism  of  a  Fatigue  Crack  Nucleation  Mechanism,"  J. 
Mech.  Phys.  Solids,  Vol.  17,  pp.  511-523,  1969. 

11.  Lin,  T.H.  and  S.R.  Lin,  "Micromechanics  Theory  of  Fatigue  Crack  Initiation  Applied 
to  Time-Dependent  Fatigue",  ASTM  STP  675,  pp.  707-728,  1979. 

12.  Lin,  T.H.  and  S.R.  Lin,  "Fatigue  Crack  Initiation  with  Creep"  Defects  and  Fracture 
Edited  by  G-C.  Shi  and  H.  Zorski  published  by  Martinas  Nijhoff  Publishers,  pp.  3-13, 
1982. 

13.  Lin,  T.H.,  S.R.  Lin  and  X.Q.  Wu,  "Micromechanics  of  an  Extrusion  in  High-Cycle 
Fatigue"  UCLA  Engineering  Rept.  88-09,  1988  and  to  be  published  in  Philosophical 
Magazine. 


17 


14.  McEviley,  J.,  Jr.  and  E.S.  Machlin,  "Critical  Experiments  on  the  Nature  of  Fatigue  in 
Crystalline  Materials,"  Proc.  Int.  Conf.  on  the  Atomic  Mechanism  of  Fracture,  Tech¬ 
nology  Press,  MIT  and  John  Wiley  &  Sons,  1959. 

15.  Mughrabi,  H,,  R.  Wang,  K.  Differt  and  U.  Essmann,  "Fatigue  Crack  Initiation  by 
Cyclic  Slip  Irreversibilities  in  High-Cyclic  Fatigue",  Fatigue  Mechanics,  ASTM  STP 
811,  pp.  5-45,  1983. 

16.  Mott,  N.F.,  "Origin  of  fatigue  Cracks",  ACTA  Metallurigica,  Vol.  6,  p.  195,  1958. 

17.  Parker,  E.R.,  "Theories  of  Fatigue",  Mechanical  Behavior  of  Materials  in  Elevated 
Temperatures,  Edited  by  J.E.  Dom,  McGraw-Hill,  pp.  129-148,  1961. 

18.  Taylor,  G.I.,  "Plastic  Deformation  of  Metals”,  J.  Inst.  Metals,  Vol.  62,  p.  307,  1938. 

19.  Thompson,  N.,  N.J.  Wadsworth  and  N.  Louat,  "The  Origin  of  Fatigue  Fracture  in 
Copper,”  Phil.  Mag.,  Vol.  1,  p.  113,  1955. 

20.  Thompson,  N.,  "Proc.,  Int.  Conf.  on  the  Atomic  Mechanism  of  Fracture,  Technology 
Press  MIT  and  John  Wiley  &  Sons,  1959. 

21.  Wood,  W.A.,  Mechanism  of  Fatigue,"  Fatigue  in  Aircraft  Structures  (Edited  by  A.M., 
Freundental),  Academic  Press,  New  York,  1956,  pp.  1-19. 


18 


Figure  3.  Creep  Strain  Distribution  of  Secondary  Slip 
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Loading:  345  c.p.m, 
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Figure  7.  Creep  Strain  Buildup  at  Free  Surface  versus  Loading  Time  N/f 


